Preconditioned methods for solving the incompressible and low speed compressible equations

Acceleration methods are presented for solving the steady state incompressible equations. These systems are preconditioned by introducing artificial time derivatives which allow for a faster convergence to the steady state. The compressible equations in conservation form with slow flow are also considered. Two arbitrary functions, alpha and beta, are introduced in the general preconditioning. An analysis of this system is presented and an optimal value for beta is determined given a constant, alpha. It is further shown that the resultant incompressible equations form a symmetric hyperbolic system and so are well posed. Several generalizations to the compressible equations are presented which generalize previous results

Fast solutions to the steady state compressible and incompressible fluid dynamic equations

For low speed flows the use of the compressible fluid dynamic equations is inefficient. The use of an explicit scheme requires delta t to be bounded by 1/c. However, the physical parameters change over time scales of order 1/u which is much larger. Hence, it is not appropriate to use explicit schemes for very subsonic flows. Implicit schemes are hard to vectorize and frequently do not converge quickly for very subsonic flows. If one is only interested in the steady state then a minor change to an existing code can greatly increase the efficiency of an explicit method. Even when using an implicit method the proposed changes increase the efficiency of the scheme. The Euler equations for low speed flows will be considered first and then incompressible flows. The method is generalized to include viscous effects. Supersonic flow is accelerated by essentially decoupling the equations

Algorithms for the Euler and Navier-Stokes equations for supercomputers

The steady state Euler and Navier-Stokes equations are considered for both compressible and incompressible flow. Methods are found for accelerating the convergence to a steady state. This acceleration is based on preconditioning the system so that it is no longer time consistent. In order that the acceleration technique be scheme-independent, this preconditioning is done at the differential equation level. Applications are presented for very slow flows and also for the incompressible equations

Flux-vector splitting and Runge-Kutta methods for the Euler equations

Runge-Kutta schemes have been used as a method of solving the Euler equations exterior to an airfoil. In the past this has been coupled with central differences and an artificial vesocity in space. In this study the Runge-Kutta time-stepping scheme is coupled with an upwinded space approximation based on flux-vector splitting. Several acceleration techniques are also considered including a local time step, residual smoothing and multigrid

Central difference TVD and TVB schemes for time dependent and steady state problems

We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy

Computation of acoustic waves in a jet

A numerical treatment of acoustic waves in a jet is described. The full time dependent Euler equations are used in both linear and nonlinear formulations. The computational region of integration is artificially bounded and boundary conditions are developed to simulate outgoing waves and to enable the computational domain to be substantially restricted. Higher order methods and coordinate transformations are introduced to further reduce the number of grid points as well as to increase the efficiency of the program. Numerical results are presented for time harmonic sources as well as for sources with more complicated time dependence

Global collocation methods for approximation and the solution of partial differential equations

Polynomial interpolation methods are applied both to the approximation of functions and to the numerical solutions of hyperbolic and elliptic partial differential equations. The derivative matrix for a general sequence of the collocation points is constructed. The approximate derivative is then found by a matrix times vector multiply. The effects of several factors on the performance of these methods including the effect of different collocation points are then explored. The resolution of the schemes for both smooth functions and functions with steep gradients or discontinuities in some derivative are also studied. The accuracy when the gradients occur both near the center of the region and in the vicinity of the boundary is investigated. The importance of the aliasing limit on the resolution of the approximation is investigated in detail. Also examined is the effect of boundary treatment on the stability and accuracy of the scheme

Orbital operations study. Volume 2: Interfacing activities analyses. Part 3: Data management activity group

A summary of the findings of the data management group of the orbital operations study is presented. Element interfaces, alternate approaches, design concepts, operational procedures, functional requirements, design influences, and approach selection are described. The following interfacing activities are considered: (1) communications, (2) rendezvous, (3) stationkeeping, and (4) detached element operations

Pseudo-time algorithms for the Navier-Stokes equations

A pseudo-time method is introduced to integrate the compressible Navier-Stokes equations to a steady state. This method is a generalization of a method used by Crocco and also by Allen and Cheng. We show that for a simple heat equation that this is just a renormalization of the time. For a convection-diffusion equation the renormalization is dependent only on the viscous terms. We implement the method for the Navier-Stokes equations using a Runge-Kutta type algorithm. This permits the time step to be chosen based on the inviscid model only. We also discuss the use of residual smoothing when viscous terms are present

Simulation of the fluctuating field of a forced jet

The fluctuating field of a jet excited by transient mass injection is simulated numerically. The model is developed by expanding the state vector as a mean state plus a fluctuating state. Nonlinear terms are not neglected and the effect of nonlinearity is studied. The results show a significant spectral broadening in the flow field due to the nonlinearity. In addition, large scale structures are broken down into smaller scales